generate total_cases_per_person =  total_cases_per_million/1000000
generate total_deaths_per_thousand = total_deaths_per_million/1000
generate total_vaccinations_per_person = total_vaccinations_per_hundred/100
generate thous_bri_per = bri_per/1000
generate unga1 =  unga*-1
generate bribery1 =  bribery*-1

# Model 1 in Table 2 #

heckprobit cn_vac unga1 thous_bri_per i.tw_ally bribery1 vaccine produce total_cases_per_person total_deaths_per_thousand total_vaccinations_per_person gdp_per_capita_thous human_development_index debt_gdp growth2020 , sel(covax_expanded = unga1 thous_bri_per i.tw_ally bribery1 vaccine produce total_cases_per_person total_deaths_per_thousand total_vaccinations_per_person gdp_per_capita_thous human_development_index debt_gdp growth2020 i.developing) vce(robust) nolog

# Model 2 in Table 2 #

heckprobit cn_vac_gift unga1 thous_bri_per i.tw_ally bribery1 vaccine produce total_cases_per_person total_deaths_per_thousand total_vaccinations_per_person gdp_per_capita_thous human_development_index debt_gdp growth2020 , sel(covax_expanded = unga1 thous_bri_per i.tw_ally bribery1 vaccine produce total_cases_per_person total_deaths_per_thousand total_vaccinations_per_person gdp_per_capita_thous human_development_index debt_gdp growth2020 i.developing) vce(robust) nolog

# Margins plots #

heckprobit cn_vac unga1 thous_bri_per i.tw_ally bribery1 vaccine produce total_cases_per_person total_deaths_per_thousand total_vaccinations_per_person gdp_per_capita_thous human_development_index debt_gdp growth2020 , sel(covax_expanded = unga1 thous_bri_per i.tw_ally bribery1 vaccine produce total_cases_per_person total_deaths_per_thousand total_vaccinations_per_person gdp_per_capita_thous human_development_index debt_gdp growth2020 i.developing) vce(robust) nolog

margins tw_ally
marginsplot, xlab(-.25 " " 0 "No" 1 "Yes" 1.25 " ", notick)

margins, at(unga1=(0(0.3)-3))
marginsplot

margins, at(thous_bri_per=(0(1)10))
marginsplot

margins, at(bribery1=(-4(0.7)3))
marginsplot